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Momentum space of a gapped quantum system is a metric space: it admits a notion of distance 
reflecting properties of its quantum ground state. By using this quantum metric, we investigate 
geometric properties of momentum space. In particular, we introduce a non-local operator which 
represents distance square in real space and show that this corresponds to the Laplacian in curved 
momentum space, and also derive its path integral representation in momentum space. The quantum 
metric itself measures the second cumulant of the position operator in real space, much like the 
Berry gauge potential measures the first cumulant or the electric polarization in real space. By 
using the non-local operator and the metric, we study some aspects of topological phases such as 
topological invariants, the cumulants and topological phase transitions. The effect of interactions 
to the momentum space geometry is also discussed. 
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I. INTRODUCTION 

Topological insulators (superconductors) have been at- 
tracting a lot of interest both theoretically and experi- 
mentally. Similar to trivial insulators, they have a band 
gap in the bulk. The striking feature that distinguishes 
topological insulators (superconductors) from trivial ones 
is that they have gapless modes on the boundaries. In 
other words, boundaries of topological insulators arc to- 
tally metallic. These gapless modes are robust under per- 
turbations and cannot be gapped without going through 
a quantum phase transition, i.e., the gapless modes are 
topologically protected. One of the best-known exam- 
ples of the topological phases is the integer quantum hall 
effect (IQHE) in 2-t-l dimensions where the Hall conduc- 
tance axy is quantized. 

Recent excitement on topological phases comes from 
the discoveries of new topological phases; the quantum 
spin Hall effect (QSHE) in two space dimensions and the 
Z2 topological insulator in three space dimensions pl-pl . 
The existence of these new phases suggests that the topo- 
logical phases are relatively generic for electron systems: 
they are not restricted to two space dimensions nor they 
do not necessarily require time reversal symmetry break- 
ing. 

Following these discoveries, the exhaustive classifica- 
tion of topological insulators and superconductors of non- 
interacting fermionic systems was proposed in Refs. p- 
pd| . The classification is based on discrete symmetries 
that are generic in quantum systems; time reversal sym- 
metry, charge conjugation and chiral (sublattice) symme- 
try. All the topological phases are characterized by cer- 
tain types of K-theory charges; Z Chern/winding num- 
bers and Z2 numbers. (For an introduction of K-theory 
charge for fermi surfaces, see [Q.) These quantities 
are defined in momentum space: the Chern number is 
associated with the Berry curvature (gauge potential), 
the winding number is given by a Wess-Zumino-Witten 
type term constructed from a projection operator, and Z2 



charge is given by some gauge invariant product of time 
reversal polarization at time reversal invariant points. 
The well-known example is again the IQHE: the Hall 
conductance axy is given by the first Chern number of 
the momentum space Berry curvature via the TKNN 
formula |13, H. For time reversal symmetric systems. 
Berry curvature docs not give a non-trivial Chern num- 
ber. Instead, the Z2 number, which basically counts the 
number of Kramers pairs or Dirac points, gives the topo- 
logical charge |iU,|,|l|. 

It is natural to ask if there are other functions or quan- 
tities in momentum space that may also play a role. It 
is known that the Berry gauge potential is not the only 
quantity that we can define in momentum space. For in- 
stance, we can generalize the idea of the Berry gauge po- 
tential to higher order tensors such as a metric. Given the 
fact that topological invariants constructed out of metric 
play a role in general relativity, it would be interesting to 
investigate physical properties of the momentum space 
metric to see if it captures new aspects of topological 
phases. The momentum space metric (the quantum met- 
ric, or the Bures metric) has considered in Refs. M, U% 
in the context of localization and polarization of Wan- 
nier functions and a relation to the cumulants. In this 
paper, we explore a connection between the momentum 
space metric and physically observable quantities. Espe- 
cially we propose a non-local operator that corresponds 
to the Laplacian in curved momentum space. We will 
also study some representatives of topological insulators 
from the point of view of the momentum space metric. 

The organization of this paper is as follows. In sec- 
tion O, wc review the construction of the Berry gauge 
potential and the Bures metric (the quantum metric). 
In section HI, we first review the known properties of 



the Bures metric and then define and investigate a non- 
local operator corresponding to the Laplacian in curved 
momentum space. In section |^ we study several topo- 
logical phases by using the Bures metric, especially for 
symmetry classes A (unitary) and AH (symplectic, spin- 



orbit). 



II. MOMENTUM SPACE GAUGE FIELD AND 
GRAVITY 

A. setup 

In this section, we will explore momentum space ge- 
ometry and topology by studying gauge field (the Berry 
gauge field) and metric (the Bures metric). Consider a 
tight-binding Hamiltonian, 

if = ^V't(r)H(r,r')V(r'), (1) 

r,r' 

where ip{r) is an A^/-component fermion annihilation op- 
erator, and index r labels a site on a d-dimensional lattice 
(the internal indices are suppressed). Each block in the 
single particle Hamiltonian H{r, r') is an Nf x Nf matrix, 
satisfying the hermiticity condition 'H^'(r',r) = H(r,r'). 
and we assume the total size of the single particle Hamil- 
tonian is NfV X NfV, where V is the total number of 
lattice sites. The components in V'(r) can describe, e.g., 
orbitals or spin degrees of freedom, as well as different 
sites within a crystal unit cell centered at r. 

Provided the system has translational symmetry. 



•H(r,r')--H(r-r') 



(2) 



with periodic boundary conditions in each spatial direc- 
tion (i.e., the system is defined on a torus T'^), we can 
perform the Fourier transformation and obtain in mo- 
mentum space 



H 



Y, ^^(k)H(k)^(k), 

kGBZ 



(3) 



where the crystal momentum k runs over the first Bril- 
louin zone (BZ), and the Fourier component of the 
fermion operator and the Hamiltonian arc given by 

V'(r) = V^ J2 e'^'^W, H(k) =^e-*-H(r), (4) 

kSBZ r 

where V is the total number of sites, and 'H(k) is an 
Nf X Nf matrix. 

The Bloch Hamiltonian 'H(k) is diagonalized by 

n{k)\u-{k))^E-{k)\u-{k)), a^l,...,Nf, (5) 

where and |w''(k)) is the a-th Bloch wavefunction with 
energy .^"(k). The fermion operator can be expanded in 
terms of the eigen functions as 



'0»(k) == ii,''(k)xa(k), i,a = l,. 



■^Nf, 



(6) 



where Xa (k) represents a fermionic operator in the eigen 
basis and given by 



X6(k)=K(k)]V,(k). 



(7) 



We assume that there is a finite gap at the Fermi level, 
and therefore we obtain a unique ground state |$(k)) for 
each k by filling all states below the Fermi level. (In 
this paper, we always adjust E'^{k) in such a way that 
the Fermi level is at zero energy.) We assume there are 
N- {N+) occupied (unoccupied) Bloch wavefunctions for 
each k with N+ + N^ = Nf. We caU the set of filled 
Bloch wavefunctions {\u^ (k))}, where hatted indices a ~ 
1,...,N^ labels the occupied bands only. The many- 
body ground state \^) is the filled Fermi sea 



JV_ 



m = 



i$(k)), i$(k))=n^i(k)io>- 



(8) 



keBZ 



a=l 



B. Berry connection and Bures metric (quantum 
metric) 

The momentum space gauge field and metric can be 
introduced by comparing two states |$(ki)) and |$(k2)) 
at different momenta ki_2- From the overlap ($(k)|$(k+ 
dk)), we define the Berry connection |1^, |l9[ A{k) = 
Afj^(k)dk^ by the expansion 



($(k)|$(k + dk)) = 1 + A^,{k)dk, 



(9) 



where ^ = 1, ■ ■ ■ ,d. Similarly, to define the momentum 
space metric, we consider the quantum distance -Dki.ka 
between two states |$(ki)) and |<I'(k2)) by I?^ ^ = 
1 - |($(ki)|$(k2))P !§. When these two states are 
infinitesimally close to each other, ki = k, k2 = k + dk, 
the quantum distance can be expanded as 



|($(k)|$(k + dk)>|2 = I - g^,,{k)dk^dk. 



(10) 



This defines the quantum metric g^^{k). 

For the case when |$(k)) is given by a Slater determi- 
nant, the overlap ($(ki)|$(k2)) is computed as a deter- 
minant of an N^ X N^ matrix made of Bloch states, 



($(ki)|$(k2)) =det5(ki,k2), 
5,g(ki,k2):=(wr(ki)|ur(k2)), 



(11) 



where |w^(k)) (a == 1, • • • , N^), are the Bloch wavefunc- 
tion for filled bands. Hence the Berry connection and the 
Bures metric are given by 



A^{k) = {u7{k)\d^u7{k)), 
5A<>^(k) = :^[{d,,u7\d^u7) + {d^uT\df,u7] 
+2{u7\df,u7){u7\d^u7] 



(12) 



(13) 



= Re (9^u. la^u- ) - (a^M- |u^ )(u. la^ii. ), 

where repeated indices are summed, and df_, = dk^, , etc. 

For the case with more than one occupied bands, one 
can also introduce the non-Abelian Berry connection by 

^"^(k) = Af{k)dk, = {u7{k)\du7ik)), (14) 



where Af,^ — — (A?,^)*. A non-Abelian generalization of 
the metric is discussed in Refs. [ pi[ . 

The above formula for the gauge field and metric can 
compactly summarized in terms the spectral projector 
onto the filled Bloch states and the "Q-matrix", which 
are defined by 

P(k) := Y. I"a (k))(«a (k)l, Q(k) := 1 - 2P(k). (15) 



Then, the U(l) part of the Berry curvature the metric 
are given by 

1 



F^,(k) = F;°(k) = -tr [Qd.Qd^Q] 



(16) 



I-G-i 5/^1^ a-nd F^^^ are the symmetric and antisymmet- 
ric part of tr [PdiyPdf^P], respectively, 2^^^ + F^^ = 
-tr [Pd,Pdf,P]. 



III. NON-LOCAL OPERATORS AND 
MOMENTS OF THE POSITION OPERATOR 

In this section, we discuss the expectation values of the 
position operator, and its second cumulant in a quantum 
ground state. For band insulators, the Berry connection 
and the Bures metric naturally enter in evaluating these 
expectation values. This is to some extent expected, since 
the definition of the Berry connection and the Bures met- 
ric involves a derivative with respect to momentum, this 
in turn suggests that it is related to the position oper- 
ator in real space. We will first review these ideas de- 
veloped in the theory of macroscopic polarization, and 
in the context of maximally localized Wannier functions 
in solids. We will then introduce yet another formalism 
for the second cumulant of the position operator, which 
is for the geometric, rather than arithmetic, mean of the 
second cumulant. We will show how it is related to the 
momentum space spectral geometry. 



A. Marzari-Vanderbilt theory 

Let us consider the following non-local operator: 



exp 



y^ia ■ rp(r) 



(17) 



where r = (xi,X2,--- ,Xd) labels sites on the d- 
dimensional hyper cubic lattice with TV* (= V) sites, and 
can be written, in terms of the lattice constant a as 



I,--- ,A^, L^aN. 
p{y) is the density operator at r, 

Nf 



p(r) = ^Vl(r)'/'i(r)- 



(18) 



(19) 



We now choose 



a = 27rn^/L, 



(20) 



where n^ is the unit vector in the ^-direction and con- 
sider the expectation value of the non-local operator. 



3m - (*k2.A./L I*) 



(21) 



The logarithm of 3 defines the expectation of the position 
operator in periodic systems 



L 

2^ 



Imln^^. 



(22) 



(The "single point formula" for the macroscopic polar- 
ization by Resta |||-|2|l.) 

For the case of Slater determinant, this can be evalu- 
ated as 



L 

2^ 



^trlnS'(k,k4-k^). 



(23) 



where k^ = 27rn^/L and the trace is taken over occupied 
states. In the continuum limit, this reduces to the Wilson 
loop of the U(l) part of the Berry connection, {x^) ^ 
In exp ^g2 tr -^(k) . 

In general, we can consider the generating function 

C(a) = (*|z„|*). (24) 

For slater determinants, this can be evaluated as 



In C{a) = ^ tr In S'(k, k + a). 



(25) 



Taking the derivative with respect to a, we generate all 
cumulants of the position operators. 



[Xi X2 



da'l^dc 



■day 



InC(a) 



Q=0 



(26) 



where the derivative with respect to a^ should be prop- 
erly discretized [M. E.g., 



(a^P^^p) 



fl'^fl/C 



a 



2ttJ NfL 



"^ ln]^dct|S'(k,k + k^)|^ (27) 



where fi is not summed. In the continuum, the second 
cumulant is given by 



\X fj,X jy J f^ 



Nfi2n) 



BZ 



dkgf,^(k) 



N 



■^i,fiv (28) 



/ 



1=1 



with V the cell volume. 

This expression for the second cumulant should 
be compared with fi = ^^ [{^'^)a — (r)a] 1 where 
(r2)a = [v/{2Tr)'']{d^Uaik)\d^Ua(M)) and (r)a = 
i\v/{2nY] J(ik(ua(k)|Vk|ua(k)), which was also used in 
p6[ . In that paper, it was showed that in a band model, 
ri is not gauge invariant and hence cannot directly be re- 
lated to an observable, while Vlj is gauge invariant. This 
gauge dependence argument would not be relevant in this 
paper. We mostly consider general continuum models. 



B. non-local operator 

We introduce yet another non-local operator rj by 



(29) 



where p is the average density, and 9(r) is assumed to 
be written in terms of a logarithm of some function, 



e(r) :=hi6'(r) 
For our purpose, we will choose 

d 



(30) 



0(r) = 
with 



-2E 



cos- 



27rnp 



2d 



a = 1/N, = l/{fV) 



(31) 



where Ne is the total number of electrons, and / = Ne/V 
is the filling fraction. By expanding the cosine in the 
decompactifying limit, Xfj,/L — > 0, 



ln6l(r)-ln^ 



^2nxn\ . r 
' = in 



and noting 



M=l 



d^r 



aN 



N 



(32) 



In- 



In 



[p(r) - P] 



= / d'^r In(rVa^) [p(r) - p] 



(33) 



where we have used the charge conservation J d'^r[p{r) ~ 
p] = 0, r] may be interpreted, in the continuum limit, as 



f f d'^r 
rj^cxpla -^ (In r^/a^) [p(r) - p] 



(34) 



In the first quantization, the density operator is given by 



p(r)=^5('^)(r-r, 



(35) 



i=i 



where r^ is the position operator of the i-th particle. 
Then, 



cxpa / d''r(ln r^)p{r) = e^^-i^"''' ^ J] (rf)".(36) 
•' i=\ 

Thus, the expectation value of 77 is nothing but the geo- 
metric mean of the position operator over TVg fermions. 
We note that by choosing 6(r) in Eq. (P9| ) properly, r] 
can represent several different non-local operators famil- 
iar in condensed matter context. When we choose 



(r) = e' «" , 6(r) = i 






(37) 



r\ is the twist operator 3 used in the theory of electric 
polarization. In two dimensions, and when we choose 



Q{v) = [z-z^r\ e(r)=ln(z-zo)' 



(38) 



where z — x + iy, and zq is some reference point, then 
77 i s the disorder operator studied by Shindou ct a l. | 25| , 
P6| . Yet another non-local operator in the form (^91) is 
discussed in Ref. [^ in the context of the entanglement 
entropy. 



C. the non-local operator for band insulators 

We now specialize to band insulators and evaluate 
the expectation value of the non-local operator 77. Be- 
low, we will evaluate the expectation value ($|77|$) 
for a band insulator. Since we can factorize 77 = 
exp U2r 0(i")p(i")] exp [— J2r ®(i")p] ' ^^6 focus on the first 
factor (will be denoted by rj for simplicity). To this end, 
we first consider 

V--C C:=exp ^\n(j){r)p{r) 

_ r 

0(r) =: [0(r)]" . (39) 

Note that in real space, 

evl(r) = 0(r)7^t(j.)^ (40) 

In the momentum basis, 

^Vl(k) = VV~' ^ 7Al(r)0(r)e+-''e (41) 

r 

For our choice of ^(r) and 0(r), 

^TZ-^k) = [ + 2d7/.|(k) 

-^(7Ai(k+k^) + 0i(k-k^))]e 

=:5]7Al(k')r(k',k), (42) 

k 

where k^ = (27r)n^/L. The matrix r(k',k) can be 
viewed as a discretized Laplacian (tight-binding Hamil- 
tonian) on the dual lattice. 

In the basis which diagonalizes H(k), 

exl(k)r' = E4(k')T(k',k)(7ir(k')|7.r(k)) 

k 

=^E4(k'ma(k',k). (43) 

k 

The matrix f'.g(k,k') (or f'!-(k,k')) is the properly de- 
fined "Laplacian" in the "curved" momentum space. 



The expectation value of the non-local operator is then 
given by 



($|e|$)-Det [f,j,(k,k'] 
($|r?|$)=Det fe(k,k': 



(44) 



where Det is taken over the hatted indices as well as the 
momentum. We have thus related the geometric mean of 
the position operator squared to the spectral geometry 
defined for the Laplacian in the curved momentum space. 



1. path-integral representation 

For simplicity, we assume there is only one band occu- 
pied. In this case, 

f,^(k,k'):=r(k,k')(ur(k)|u7(k')> 



2dSu> 



k',k 



a 



k+ky^M 



U,{^)+K'm-LUI{^) 



k',k-k,,'^A'"> 



^L=l 



2dl - X, 



(45) 



where we have assumed two bands a = ± 1 , and we fill the 
lower band, a ~ —. The link variable U^{\i.) = (u(k + 
kp)|w(k)) is decomposed into the amplitude (hopping) 
and phase (gauge field) parts as 



t^{k) = |(u(k + k^)|u(k))|, 

gA,(k) ^ (»(k + kM)l"(k)) 
\{u{k + k^)\u{k))\' 

For an infinitesimal k^. they are approximated as 

t^(k) ~ exp f --5^^(k)Afc^AA:^ j , 

e^.(k)^exp(^^(k)AfcJ 



(46) 



(47) 



where /i is not summed, 

In order to compute DetT = Det (2dl — X), we con- 
sider the following formal expansion in terms of y = 
X/{2d): 



Det (/ - Y) = exp [Tr ln(/ - Y)] 



exp 



CXD _, 



(48) 



For a given n, the trace is represented as a sum over 
paths (loops): 



Tr {Y''' 



\Cp\=n 

EE n u,ik) 

p L-'p Op 

|Cp|=n 

EE n *.w 

P c„ c„ 



,^f.(k) 



(49) 



where p G BZ is an initial point for paths, Cp represents 
a path starting from p and going back to p with n steps, 
and He ^t^i^) i^ ^^^ path-ordered product. Thus, the 
expansion can be interpreted as a summation over closed 
loops. It should be noted that He ^t^i^) i^ actually does 
not depends on the initial point p at which one starts to 
draw a path. 
Thus, 



1 / 1 



|Cp| = 



ln($|C|$)=lnDetf^^- - Y. E ^"'- 

n=l ^ '^ P Cp 

where S=- g^,^,{k)Ak^,Ak^, - A^Ak^,.{50) 

«/ O p ''Op 

We can see that for a smaller value of the metric g^,y, 
the non-local operator takes a smaller value. Especially 
the trivial metric g^^, ~ gives the minimum expectation 
value of the operator. 



IV. BURES METRIC FOR TOPOLOGICAL 
INSULATORS AND SUPERCONDUCTORS 

In this section, we explore momentum space geome- 
tries and topologies of gapped phases, applying gen- 
eral formalism presented in the previous sections. In 
particular, we will consider examples of topological in- 
sulators/superconductors constructed in Ref. |lfl|. In 
Rcf. jl^, gapped Hamiltonians of Dirac type are sys- 
tematically constructed, which realize topological insula- 
tors/superconductors in all symmetry classes and in all 
dimensions. For these, we will discuss the momentum 
space metric and the second cumulant of the position 
operators, by following the Marzari-Vanderbilt formal- 
ism, and also by following the formalism developed in 
Sec. 



IIIB 



We will also discuss topological invariants made of the 
Bures metric, such as the Euler class (d = 2n), the Pon- 
trjagin class (d = 4n), the gravitational Chern-Simons 
invariant (d = 2n -|- 1), and the dimensionally continued 
Euler density {d = 2n + l). The dimensionally continued 
Euler densities have the following forms 



r^^i 



AE^'■' 



'^tJ-l---P2n 



(51) 



where rj'^i--^^^ = f^MiM2 ^ . . . ^M2„-iM2„ jg ^he wedge 
product of the curvature two form and ii^^sn+i-Md ^ 
j^fi2n+i /\ . . . iJA'd is the wedge product of the co-frames of 
the orthonormal frames. Among the dimensionally con- 
tinued Euler densities, there is one which is a topological 
invariant for a given dimension d. 

Following Ref. [[L0| , we start with topological Dirac in- 
sulators (Sec. IV A) in d = 2n- and d = (2?! + 1)- dimen- 



sions and discuss the general properties of the momen- 
tum space metric. These Dirac representatives realize 
topological insulators/superconductors characterized by 
an integer topological invariant for all symmetry classes 



and for all dimensions (sometimes called "primary se- 
ries"). We will then focus on specific examples of the 
Dirac representatives in d = 2,3,4, which include the 
QHE, the QSHE, and the 3D Z2 topological insulator. 



A. Dirac representatives: generalities 

1. odd space-time dimensions, d = 2n 

We start with the d = 2n-dimcnsional topological 
Dirac insulator which is defined in momentum space by 



d=2n 



'^(2n+l)(^) — Z^ ^a'rt2n+l) +''^'^{2n+l)- (^'^) 



a=l 



Here, ka=i ...,d ^^^ a component of the d-dimensional mo- 



mentum, and r' 



a=l,...,2n+l 
(2n+l) 



are (2" x 2")-dimensional her- 



mitian matrices and satisfy {r?2n+i)' ^(2n+i)} ^ 2(5a,6. 
The two eigenvalues of the Dirac Hamiltonian are 



±X{k), \{k)^ Vk^ + m^, 



(53) 



where k = |k|. 

The Dirac Hamiltonian ( 52 ) represents the topological 
insulator/superconductor in class A (d = 2n), which is 
characterized by lack of any discrete symmetry, and by 
an integer topological invariant. For a given d = 2n, the 
Hamiltonian ( ^2[) respects either one of TRS or PHS. For 
this reason, the Dirac Hamiltonian (^2|) also represents 
the topological insulator/superconductor in class AI (d = 
8m), class D (d = 8m + 2), class AH (d ~ 8m + 4), class 
C (d = 8m + 6), where m G I^. 

The momentum space metric for the Dirac Hamilto- 
nian (^2|) is given, in the Cartesian and in the radial 
coordinates, by 



ds' 



-^ {Sab>^^ - kakb) dkadkb 
,2 p 



A4 



■dk^ 



X 



:dQd-i 



(54) 



respectively. Here, ds^ is the line element in momen- 
tum space, and dQd-i is the line element of the (d — 1)- 
dimensional unit sphere. 

Near the origin fc = 0, and near fc = 00, the metric is 
given by 



ds' 



2" 



^{dk^ + k^dfld-i), fc^O 



m^-— + dfld-i ) , k —)■ CO. 



(55) 



The metric for fc — >■ is equivalent to that of a sphere, 
while the metric near fc = 00 can be brought into a form, 
with z = 1/fc, ds^ ~ ^{m'^dz^ +d^d-i)- I-e., the met- 
ric of a flat geometry. Therefore the momentum space is 



smoothly capped at fc = while it has a finite radius at 
fc = 00 (2: = 0). This shows that the geometry of momen- 
tum space, defined by the Bures metric, is a semisphere 
or a cigar. 

a. regularization For the Dirac Hamiltonian (p2), 
the behavior of the wavefunctions at high energies is non- 
trivial. This is signaled by the fact that the Chern-integer 
for d = 2n, when computed for the Dirac Hamiltonian 
(p2|), is not quantized in integer unit but in half-odd in- 
teger unit. In turn, this is closely related to the fact 
that the momentum space geometry detected by the Bu- 
res metric is that of a semisphere. For these reasons, 
while the Dirac Hamiltonian ( p2| ) is capable of describ- 
ing a transition between two phases with different topo- 
logical charge (as we tune m), it does not represent, in 
a well-defined fashion, a topological or a non-topological 
insulator by itself |10[| . 

If we wish, such non-trivial behavior of the wavefunc- 
tion can be regularized, by modifying the Bloch Hamil- 
tonian at large momentum, for example, by replacing the 
constant mass by 



m{k) ~ m ~ Ck 



(56) 



where C is some constant. With the regularization, the 
71-th Chern integer is quantized in integer unit, Ch„ = 
1, for sgnm = zbsgnC, respectively (d = 2n). 

The metric for the regularized Dirac Hamiltonian is 
given by 



ds^ = — 



—^(ifhadf^a T 

{m + Ck^f 
A4 



(4mC- 1) 
A4 



kaf^b^'^a^^b 






dk^ + '—dVld^i 



(57) 



where A(fc) = •\/fc^ + m^(fc). 

Near the origin fc = 0, the metric shows the same 
asymptotic behavior, ds^ ~ ■§^{dk'^ + k^dVtd-i), as in 
the unregularized case, while near fc = 00, 

fdk^ 1 jr. ^ 



ds^ 



2" 
8C2' 

2" 
8C2' 



fc4 fc2 



.^{dz'^ ^z'^dUd^ 



(58) 



Therefore momentum space is smoothly capped both at 
fc = and at fc = 00 (z = 0). This shows that the 
topology of momentum space is a sphere. 

It is also interesting to see the difference between the 
regularized and unregularized models at criticality. In 
the limit m — > 0, the metric for the unregularized model 
( p4| ) reduces to ds'^ = ^dVid-i, i.e, g^k = 0. On the 



other hand, for the regularized model, (57), the metric 
reduces 



ds^ 



2" 



C^dP 



da 



d-l 



(l-hC2fc2) 



(59) 



l-HC2fc2 

1 does not vanish at 



Notice that the coefficient of dild 

k — 0. Therefore, the topology of momentum space 

changes to create a new "hole" at the origin. 



Alternatively, one can put the theory on a lattice to 
regularize. This amounts to consider the following lattice 
version of the Dirac Hamiltonian: 



d=2n 



W(2=n+i)(k) = E sinfcar^2„+i) + /(k)r^^'++\) . (60) 



a=l 



where ka G (— tt, +7r], and /(k) = {m + d) — 
2C^^~^"cosfca. The eigenvalues are given by ±A(k) 
with 



Xik) = Jj2 (sinfca)2 + /2(k). 



(61) 



near k ^ 0, one recovers the Dirac 



By expanding 

Hamiltonian in the continuum, (52) and (|56|). The lat- 
tice Dirac Hamiltonian (pG) will be used to discuss the 
expectation value of the non-local operator rj. The metric 
in the Cartesian coordinates is given by 



gabO^) = -^-r^ [cos ka cos kbSab + 4C^ sin ka sin kb] 



2" 



sin kr, sin kh 



X [cosfca + 2C/(k)] [cos kb + 2C/(k)] , (62) 

where repeated indices are not summed over. Observe 
that the volume element ^/g can be zero, say, at ka = 
±7r/2. 



2. even space-time dimensions, d = 2n — 1 

A Dirac Hamiltonian in even space-time dimensions 
(d = 2n — 1) can be obtained from ( p2|) by formally re- 
placing one component of the momentum by a mass term, 
^211 — > rni (dimensional reduction), 

d=2n~l 



^(2n+l) (^) — Z^ ^^r^jn+i) 



a=l 



+rnir^2Vi) + ™2r^r„+i) 



(63) 



When one of the masses is set to zero, the Hamiltonian 
anticommutcs with ^(2n+i) i ''^^'^ ^^^^ i^ ^ topological 
insulator in class AHI, i.e., there is chiral (sublatticc) 
symmetry In addition to this chiral symmetry, for a given 
d = 277. + 1, the Hamiltonian (pSh respects either one of 
TRS or PHS. For this reason, the Dirac Hamiltonian (|6^ ) 
also represents the topological insulator/superconductor 
in class BDI (d = Srn +1), class DHI {d = 8m + 3), class 
CH (d = 8m -I- 5), class CI (d = 8m -I- 7), where m g Z. 

With the Dirac representatives (^2|) and (|63| ) , we cover 
all topological insulators/superconductors labeled by an 
integer topological invariant (called "primary series"). 
Topological insulators/superconductors characterized by 
an binary topological invariant (Z2) or an even inte- 
gral topological invariant (2Z) can be derived from these 



by dimensional reduction. We will not give an exhaus- 
tive study of the momentum space metric for these "de- 
scendents" (i.e., topological insulators/superconductors 
characterized either by Z2 or 2Z topological invariant). 
Rather, we will focus on the QSHE (class AH in d = 2) 
and the 3D Z2 topological insulator (class AH in d = 3) 
(see below). 



The two eigenvalues of the Dirac Hamiltonian (63) is 

(64) 



±A(fc), A(fc) = y^fc2 ^ ,„2 _^ , 



The calculations of the metric, etc., go parallel with the 
case oi d = 2n. The quantum metric is given by 



gabCk) = —J {SabX^ - kah) 



(65) 



This is the same metric as that of the d = 2n case. There- 
fore the geometry of momentum space is a cigar. 

As before, non-trivial behaviors of the Bloch wavefunc- 
tion can be regularized, by replacing the constant mass 
by 

mi^2 ^> "^l,2(fc) = m-1,2 - C'i,2^^- (66) 

With the regularization, the metric can be written 

— '^y^dkadka + , 2 " i> 



ds^ 



4|q2|m| 



(2m • C - 1 ) 



dkadkb 



oA 



8A2 



k^dn 



d-l 



(67) 



where A(fc) = V^^ + m ■ m, and we have introduced a 
vectorial notation, for instance, C ■ m ^ Cimi + C2m,2, 
|C|2 = C-C, etc. 

The structures of the metric at fc = and k — 00 are 
the same as those of d = 2n case. Therefore, the topology 
of the momentum space is a sphere. 



B. example in d = 2: the IQHE (class A) 



Consider three mutually anticommuting, hermitian 



matrices V 



,a=l,2,3 

(3) 



{cx, Oy, CTz}, where o^^y^z are the 
2x2 Pauli matrices. We consider the following (unregu- 
larized) d = 2-dimensional massive Dirac Hamiltonian, 



'Hl)^(k) = k^cFx + kyffy + maz- 



(68) 



This is a class A Hamiltonian. 

Assuming m > 0, the Berry connection and the Berry 
curvature are given by 



A,(k) 



Z/C-i/ 



2A (A -F m) ' 
F^y (k) = dk^Ay - dk^A^ 



A,(k) = - 






2A (A + m) '' 



im 
'2A3' 



(69) 



The Chcrn number is non-zero 

i 



Chi 



. d K> i^ xy 



1 m 

2 M 



(70) 



which is nothing but the Hall conductivity axy This 
model can be realized in the low energy limit of the hon- 
eycomb lattice model introduced by Haldane [pSl . While 
the Chcrn integer in the unregularized model (pq) in 
quantized in half-integer unit, with regularization, it is 
quantized in integer unit. 

a. the Marzari- Vanderbilt cumulant We now discuss 
quantities constructed from the momentum space met- 
ric. Let us first discuss the Marzari- Vanderbilt cumulant 
^i,fj,u- Plugging the concrete form of the metric, we ob- 
tain, for the diagonal component. 



^I,xx — ^i 

1 

~ 16^ 



yy 



(71) 



(2Cm - l)fc2 - 2m2 
2(C2/i;4 _ (2Cto - l)k^ + m2) 



2Cm-|-l /2C2fc2-2Cm-hl 

: arctan 



— s- 



y/ACm - 1 
1 

1 2Cm + 1 



32 V4Cm - 1 



V \/4Cm - 1 

/-2Cm+l 
WiCm - 1 



/C— CXD 



fc=0 



1 — arctan 



while the off-diagonal component is vanishing fli_xy = 0. 
The behavior fli_xx as a function of the mass param- 
eter m is shown in Fig. Il|. Observe that the plot is not 
symmetric about m = 0: the cumulant for the positive 
mass m > is larger than for the negative mass tti < 0. 
This shows that the electrons tend to be delocalized in 
the topological phase (m > in this case). In particular, 
the cumulant goes to zero as m — >■ — cxd, while it increases 
for positive m. This behavior should be compared with 
the lattice calculation in Refs. E^, p9], where the cumu- 
lant remains finite in the topological phase. Had we not 
regularized the Dirac model (pSh, fli^xx as a function of 
the mass parameter m would be symmetric with respect 
to m = 0. At the critical point m = 0, the cumulant 
diverges logarithmically 



logm. 



(72) 



b. topological invariant The topological number for 
a closed 2 dimensional manifold is the Euler number 



X 



1 

An 



(fk^R, 



which is equal to 2 — 2g — 6 for a smooth manifold with 
genus g and boundary number b. 

As we mentioned above, the momentum space is a 
cigar and non-compact if we do not introduce a regulator. 
With the regularization 7TT, — )- fh{k) = m—Ck^, the topol- 
ogy of the momentum space is a sphere and so the Euler 
number is x = 2. Note that when sgn(r7i) = — sgn(C), the 
metric (p7|) has a coordinate singularity at k^ = —m/C. 




FIG. 1: The Marzari- Vanderbilt cumulant Q.i,xx for the 2D 
Dirac model (pS) as a function of the mass parameter m. 
The mass term is regularized by replacing m in ( pq) by m ^> 
fh{k) = m — Ck^, and we have set C = 1. 



c. second cumulant: the geometric mean We now 
study the expectation value of the non-local operator 77 
for a topological Dirac insulator in class A. Here, we con- 
sidered the lattice regularized model ( pOJ ) with d — 2. 
The Dirac model describes the chiral p-wave supercon- 
ductor (p-wavc SC) on the square lattice defined by the 
tight-binding Hamiltonian 

H = Y, ^(r)t (^ _\ ^^ ) iiv + i) + h.c. 
t iA 



+'^('")^ i iA '" )^(r + y) + li.c. 



~fi 



V^r), 



(74) 



where the multi-component fermion annihilation oper- 
ator at site r, '0(i")j is given in terms of the electron 
annihilation operators with spin up and down, c* 1 (r) , 

= (c.^(r),c](r)), and we take t = A = 1 and 



as ifi (r 

fi — m + 2. The chiral p-wave SC has been discussed 
in the context of superconductivity in a ruthenate and 
paired states in the fractional quantum Hall effect pO[ - 
|3q| . There are four phases separated by three quan- 
tum critical points at /i = 0, ±4, which are labeled by 
the Chern number as Ch = {\^\ > 4), Ch = -1 
{-4: < fi < 0), and Ch = +1 (0 < /i < +4). The 
non-zero Chern number implies the IQHE in the spin 
transport [p2| . 

In Fig. H, we plot the expectation value of the non- 
(73) local operator, a^^ In ($|77|<I>), obtained by numerically 
diagonalizing the Laplacian operator in momentum space 
f,^(k,k'), y-iln($h|$) = V-ilnDct[f,^(k,k')]. 
The behavior of V ^ ln($|77|$) as a function of /i is qual- 
itatively similar to the Marzari- Vanderbilt cumulant: (i) 
it is larger in the topological phase (—4 < /i < 4) than 
in the non-topological phases (|/i| > 4), and (ii) near 
the critical points /i = ±4, 0, it shows a peak structure. 
Due to the system size limit accessible by numerics, it 
is admittedly difficult to see if V~^ \n{^\rj\^) develops a 



t:^ 



bJO 

o 



1.260 



1.255 




FIG. 2: The expectation value of the non-local operator 
V^^^ In ("1>|?7|$) for the lattice Dirac Hamiltonian (lattice chi- 
ral p-wave superconductor) (fTJ) as a function of /i. 



singularity at the critical points. 

Some technical comments are in order. Computing 
T^~^ ln($|77|$) by numerically diagonalization the mo- 
mentum space Laplacian, and then summing over log- 
arithm of eigenvalues can be quite tricky, since some 
eigenvalues are very close to zero. In order to obtain 
a meaningful result, we can modify (or regularize) 0(r) 
as 



= ^25:, 



d 



N 



2d 



-2^eos '"^"^ + ^"^^ +2(rf + Jrf). (75) 



^L=l 



N 



Here, adding 6d<^ d shifts all eigenvalues £i of the Lapla- 
cian, Ei — > e -|- 25d, thereby prevents having an eigen- 
value which is too close to zero. The expectation value 
In Det [T-g(k, k')] can be evaluated for small but finite 
5d^ and the results can be interpolated to (5d — > 0. Adding 
< 5n^ < 1 amounts shifting the "location" of the non- 
local operator at which it is inserted; the function 0(r) 
is centered at a site r = 0, but alternatively, one can de- 
fine ^(r) in such a way that it is centered at the center 
of a plaquette (i.e., a dual site), (5n^ = 1/2. Or, more 
generally, one can take the center of 0(r) anywhere in a 
plaquette, < ^n^ < 1. In Fig. |[ the expectation value 
In Det [T.jj(k, k')] is averaged over different centers (dif- 
ferent choice of < Sn^ < 1) to obtain smooth curves. 



C. example in d = 2: the QSHE (class All) 

Consider five mutually anticommuting, hermitian ma- 
trices 



r' 



a=l,...,5 
(5) 



{a^, ay, ttz, l3, -i/37^}. 



where we are using the Dirac representation. 



a, 
cJi. 



1 
-1 



1 

1 



(77) 



Observe that axayOizP = — i/375. We will use these 
gamma matrices to construct Dirac representative for 
d = 2,3,4 below. 

Let us consider the following d = 2-dimensional Dirac 
Hamiltonian. 



Jif^-^'^ik) = k^ai + kya2 
fh{k) = 771 — Ck . 



"75 + m{k){-il3j5 



(78) 



This Hamiltonian is time reversal invariant, 

T-H*(-k)T-i =-H(k), T = iTx(g)ay. (79) 



The Dirac Hamiltonian of type (|7^) was used in Refs. 
[|l[ ^ to discuss the QSHE. In the Kane-Mele model on 
the honeycomb lattice, the mass m is generated by a 
time-reversal spin-dependent second neighbor hopping, 
whereas the a-term arises from the Rashba spin-orbit 
coupling. Time reversal invariance is crucial in the 
QSHE. When the edge states consist of even number 
of Kramers' pairs of counter propagating modes (heli- 
cal modes), they can be gapped under perturbations. On 
the other hand, when the number of helical modes is odd, 
they cannot be gapped unless time reversal symmetry is 
broken. 

The eigenvalues of the hamiltonian (^S) are 



E{k) 



±^t^ 



±Xl 



V'(fc±a)2-H(m-Cfc2)2 



(80) 



Notice that for a = 0, there is a critical point at tti = 
and fc = 0. For a > 0, m/C > 0, there is a critical point 
at rn/C = c? and k = yJm/C = a. As observed in Refs. 
p6| , |37[ |, there is a metallic phase near m = for a 7^ 
sandwiched by insulating phases, while for a = 0, the 
metallic phase shrinks to a critical point. 
The Bures metric is given by 



ds^ 
1 

4 



{m + Ck{k-2a)f {m + Ck{k + 2a)f 



1 + 



A-{kf A+{kf 



v/A+(fc)A_(fc) 



dip\ 



where 



A±{k) := (m - Cky + {k±a) 



dk^ 
(81) 

(82) 



Observe that gkk is a non-negative function as seen from 
(pl|), while the behavior of g^^ is different for topologi- 
cally trivial and topologically non-trivial phases; g^p^p is 



(76) zero at fc = yJrnjC if 



10 




at the critical point m = Ca^{= 2 in Fig. ||), the cumu- 
lant takes small but finite value. The cumulant changes 
very rapidly around the critical point m — Cc? . Above 
this value a = 1, the cumulant monotonically increases. 
Another important property of the cumulant is that for 
fixed a, the cumulant takes larger value in m — Co? > 0. 
As studied in [g, ^ , m — Ca^ > region is the Z2 topo- 
logical superconductor phase. Therefore, this shows that 
the electron tends to be delocalizcd in the topologically 
non-trivial phase. 



FIG. 3: The Marzari-Vanderbilt cumulant i^i^xx{= O^/.y,,) for 



the two-dimensional Z2 topological Dirac insulator ([781) as a 
function of the mass parameter m. Here, C = 2, and a — 
0, 0.2, 0.4, 0.6, 0.8, 1, 1.2, 1.4, 1.6. from the left to the right. 



D. examples in d = 3: chiral topological insulator 
(class AIII) and Z2 topological insulator (class All) 

As in the general discussion before, we can consider 
two mass terms in d = 3, 



Let us first discuss several asymptotic behaviors of the 
metric. We first consider the case of m 7^ 0. Near the 
origin fc = 0, the metric is given by 



ds' 



2(m? + a^) 



2\2 



[dk^ + k^difi^) + 



(83) 



On the other hand, asymptotic value of g^k {k — > cx)) 
is nQiu ■ By making a coordinate transformation z = 
l/(-\/2CA;), the metric is given by 



ds — dz 



z'^dip 



(84) 



Both ( p3j ) and ( |84| ) are the metric of a sphere, and 
therefore we see that the momentum space manifold is 
smoothly capped both at the origin and the infinity. 

On the other hand, when m = and a 7^ 0, the metric 
near the origin fc = is 



ds^ 



C^ 
2a2 



2n2 



d(/c^) 



(fc 



2\3 



-dip' 



(85) 



This metric has a singularity at the origin. This should 
be compared with the case of m = and a = 0, for which 
the metric near the origin is 






^v. 



(86) 



The ip circle does not shrink to zero at the origin and the 
metric is analytically continued to the positive energy 
region. Note that the energy becomes zero at fc = and 
the gap is closed. 

d. the Marzari-Vanderbilt cumulant The Marzari- 
Vanderbilt second cumulant is calculated numerically in 
Fig. y. At the critical point a — m — 0, the cumulant di- 
verges. For finite a, the cumulant is finite with a peak at 
a,t m = Ca^, whose peak value increases as a increases. 
In another critical case (pink line in Fig. 0) where the 
radial direction of the momentum space shrinks to zero 
at certain value of fc, the cumulant gets two peaks and 



^(5f(k) 



z^ ^^r^g 



miF^g) 



mzF^g). 



(87) 



When cither one of the masses are switched off, this 
Hamiltonian belongs to the symmetry class AIII, or AIL 
As discussed in Ref. 0] for class All and in Rcf. [Q for 
class AIII, non-trivial topology of the Dirac Hamiltonian 
is signaled by the existence of the 0-term in the effec- 
tive action for the linear electromagnetic response with 
6 ^ TT. The non-zero (and quantized) 0-term, in turn, 
is related to the topological invariant made of the Berry 
connection, which is given for the Dirac model (pW), 



CS3 = ^ / d^fce^'^''tr(A,AA, 



1 m2 



arctan 



■K |m2| 



i?)V2 



—A^A^Ap 

-.1/2 



mi 



{m\ + mfy^^ + mi 



In particular, when mi = 0, CS3 — (l/4)m2/|?7i2|- As 
the CS term is determined only modulo 1, the fractional 
part (l/4)sgn(m2) is an intrinsic property. On a lat- 
tice, as there are two-copies of the 4x4 Dirac Hamil- 
tonians, the total Chern-Simons invariant is given by 
CS3 = (l/2)sgn(m2), and exp(27riCS3) = -1. 

Let us now set mi = 0, and, as before, regularize the 
mass, 



^(5) (k) = k^ai + kya2 + fc^aa + m{k){~if5^^), 



m{k) 



Ck^ 



(89) 



This Hamiltonian can be obtained from the Dirac repre- 
sentative of the QSHE by replacing a in the two dimen- 
sional hamiltonian (Pq) by fc^. 

The metric for the Dirac model (py) is 



ds' 



1 



2[(Cfc2-TO 

+k^{de' 



fc21 



(Cfc2-Hm)2 
(Cfc2 - m)2 + fc2 



dfc2 



sm 



(90) 
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FIG. 4: The Marzari-Vanderbilt cumulant Q.i^xx(y= ^i,yy) 
as a function of the mass parameter m for for the 3D Dirac 
model (H) with C = 2. 



where we used the spherical coordinate, {kx,ky,kz) = 
{k sin 6 cos (f>, k sin 6 sin 0, k cos 6*) . 

The behavior of the metric at fc = oo is similar to that 
of the d = 2 case. Here, we check the regularity of the 
metric at the origin. For m ^ 0, the metric is 



ds^ 



1 



2m2 



[dk^ + k^ {de'^ +sm'^ 



while for 771 = 0, C 7^ 0, it is 



(72 

ds = — dk" 
2 



-(d9' 



sm 



'■)■ 



(91) 



(92) 



The momentum space is regular for finite m, while at the 
critical point m = there is a hole at the origin fc ~ 0. 
Note that when sign(C) = — sign(TO), there is a coordi- 
nate singularity gkk = at Ck^ + m = 0. This suggests 
that the cumulant for sign(C) = — sign(m) phase would 
take smaller value than that for sign(C) = sign(m) phase. 
We will show this explicitly in the next paragraph. 

e. the Marzari- Vanderbilt cumulant The diagonal 
component of the Marzari-Vanderbilt cumulant ^i^xx{= 
^i,yy) for the 3D Dirac Hamiltonian (p9|), obtained by 
numerically integrating the metric over three dimensional 
momentum space, is shown in Fig. 0. The off-diagonal 
component is identically zero, ^i^xy = 0. The cumulant 
does not diverge at tti = 0, while its first derivative with 
respect to the mass m changes discontinuously: 



*'/,a:a; 



487rC 247r 



1 



487rC 127r 



m + 



-0 



-0 



(93) 



As mentioned above, the negative mass (tti < 0) phase 
has a coordinate singularity {gkk = 0) and the cumulant 
becomes smaller compared to the positive mass phase. 
The cumulant gets larger for large positive ttt. However, 
since this would come from the UV region of the integra- 
tion, it would not have physical significance. 



/. topological invariants In d = 3 momentum space, 
one can think of two invariants. Both of them are related 
to a Chern-Simons characteristics class, derived from the 
Pontrjagin, and the Euler characteristics classes, respec- 
tively. 

The Chern-Simons invariant derived from the Pontrja- 
gin class is the gravitational Chern-Simons term, which 
was discussed, for example, in Refs. [|o[ |3^, and is given 
by 



^GCS = ^ / '^^^^^ (^^^ + 3^^ 



(94) 



where a; is a connection one form. This invariant turns 
out to be zero, since the integrand vanishes identically 
for the metric (p(\). 

Another invariant is the dimensionally continued Euler 
number. In three-dimensions, dimensionally continued 
Euler number can be written as 



1 



2(47r) 



eabc^"" A e 



where eabc^"'' Ae^ = eabcR^'kie'' Ae' Ae^ 



(95) 



m-'abVad^k. 



For the unregularized model, the Dirac (^ with C = 0, 
the dimensionally continued Euler number is given by 



X 



3^2 



(96) 



It is straightforward to check that this result still holds 
for the regularized model where ttt,2 is replaced by 
fn2{k) ~ 1712 — C2k^- When ttti = and ttt2 = 1 with 
6*2 = 1, the integral can be evaluated analytically, and 
gives the same result as above. 



E. example in d = 4 

The gamma matrices (pO) can be used to construct a 
d = 4 dimensional Dirac Hamiltonian, 



^(sfW^E^-'^r^^ 



(5) 



™r^5). 



(97) 



The Hamiltonian (B^ is a topological insulator in class 
A with the non-zero second Chern number constructed 
with the Berry curvature [R, E^|. The Hamiltonian ( p7| ) 
respects time-reversal symmetry which squares to be -1, 
and hence belongs to symmetry class AH. The 3D Dirac 
Hamiltonian (|87|), which is a Z2 topological insulator in 
class AH, can be obtained from ( P7| ) by dimensional re- 
duction. M. 

g. topological invariants The Euler number calcu- 
lated by using the Burcs metric is 



1287r 



l^Jd'kVge 



klmn abed T) p i 



(98) 
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On the other hand, the Pontrjagin number (the gravita- 
tional instanton term, or the signature r) is identically 



^=3i/'^'^v^^" 



^lr^ 



rnnab jjlr 

^ ab 



0. (99) 



V. CONCLUDING REMARK AND 
DISCUSSION ON INTERACTION EFFECTS 

In this paper, we studied several properties of quan- 
tum ground states in momentum space by using the Bu- 
res metric. The Berry gauge potential and the Bures 
metric naturally appear in the expectation value of the 
non-local operator which measures the charge distribu- 
tions. Namely, the Berry gauge potential measures po- 
larization in real space and the Bures metric measures 
the second cumulant. These are the first two expansions 
of the generating function (p3) . 

Wc introduced the distance square operator in real 
space which corresponds to the Laplacian in the curved 
momentum space. In the path integral description, this 
operator is given by the summation of geodesic lengths 
and phases of all closed loops in momentum space. 

These operators, the second cumulant and the distance 
square operators, are applied to analyze some topologi- 
cal insulators. In certain phases, the Bures metrics pos- 
sess the coordinate singularities where some components 
of the metric vanish. As a result, the cumulants take 
smaller values in those phases compared to other topo- 
logically non-trivial phases. These operators also capture 
the singular behaviors of the quantum phase transitions. 

We also calculated several topological quantities such 
as the Euler number by using the Bures metric. These 
are analogous of the Chern number for the Berry gauge 
potential, which does characterize topological phases and 
has physical meaning. At this point, it is not very clear 
if the topological numbers for the Bures metric directly 
related to some observational quantities. The Euler num- 
ber, for instance, measures the topologies of the momen- 
tum spaces. Naively these are spheres or tori depending 
on regularizations. However, we have seen in the case of 
the 2d Z2 topological insulator that the topology could 
be non-trivial: one cycle shrinks to zero at finite momen- 
tum k. Also, we have seen that there are several types of 
singularities in the momentum space geometries for both 
of the Z and Z2 insulators. We would leave the investi- 
gation of observational consequences of these phenomena 
for the future works. 

Our analyses so far are based on free theories. It is 
interesting to sec how the interaction effects change the 
story. In interacting systems, while the ground state can- 
not be decomposed as a product the Slater determinants 
as Eq. (^), the Berry connection and the Bures metric 
can be defined, in principle, by twisting the boundary 
conditions B^. 

A yet another possible approach for interacting sys- 
tems is to note that the spectral projector Q(k) is nothing 



but the Fourier transform of the equal time correlation 
function ($|^|(r)^^(0)|$) Q. A proper generalization 
of the momentum space metric in interacting systems can 
then be defined in terms of the equal time correlation 
function. 

As an example, let us consider the (2+l)-dimensional 
Gross-Neveu model. I.e., the model with the Dirac ki- 
netic term as in (|6q ) , and with a four-fermion interaction. 
To the leading order in the large-TV expansion (here, N 
is the number of fcrmion flavors), the Euclidean fermion 
propagate is evaluated as (see, for example, Ref. Gsl) 



Gik^) 



ilfikf" 



m 



fc2 



(100) 



where /i ~ 0,1,2,3, fc^j=o,i.2,3 arc the space-time mo- 
menta, A is the ultra-violet cutoff, and a ~ 0{1/N). 

The equal time correlator can be read off from Q by 
/_ dkoQ{ko,ki) = G{ki). With the equal time corre- 
lator, we can compute the large-A'^ corrected quantum 
metric. The calculation amounts to modifying the pro- 
jector Q by Q(k) = •H(k)/A(k) ^ g(k) = •H(k)/A(k), 



where A(k) 



.2\2 



, and A(k) 



fc2 



The exponent 7 is 0{1/N). The metric is then given in 
the radial coordinate by 






1 



(7^ - l)e- 



A2+27 



A4+27 



dk^ 



PdVld- 

A2+27 



In the limit m — > 0, this becomes 



ds^ = — 



72 „ 1 



(101) 



(102) 



The momentum space is singular at fc = as in the case 
of non-interacting systems. However, the structure of 
the singularity is different. It would also be interesting 
to notice that in the absence of interaction, the metric 
is not deformed as we change the number of flavors TV. 
On the other hand, in the interacting case, the metric 
is deformed as a function of TV. This would also lead to 
deformation in the second cumulant. 

Another way of considering the interacting systems is 
to look at string theory construction of the topological 
phases. As pointed out in |l^ and constructed explicitly 
in E3, K-theory charges for the topological insulators 
(superconductors) are closely related to those of D-branes 
in string theory. Since D-brane systems are interacting 
ones by its construction, it would be interesting to con- 
sider momentum space geometries from D-brane point of 
view. 
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